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Parrinello and Rahman l] proposed a molecular dy- 
namics (PRMD) for crystal structure under external con- 
stant pressure in 1980. Many simulations of structure 
phase transitions induced by external forces have been 
carried out with it even till recently 0, E3- [J 01 - espe- 
cially after it was combined with the well-known Car- 
Parrinello MD. It uses the scaled particle position vectors 
Si T — (£i,r]i£i) as generalized coordinates, defined by 


with the net force Fj acting on particle i. In order to 
interpret the Eq. (2) of the PRMD's original paper[IJ ac- 
curately, especially the force term, let us deduce it again 
in the following three case steps. 

Case 1, from the Lagrangian dynamical equation, we 
can easily reproduce Eq. JSJ) by using as "general- 
ized coordinates" and using the complete kinetic energy 
Eq. to replace K' in Eq. (gj, with 


= £ja+77ib+CiC = Hs, 


(1) 


where is the ith particle position vector, a, b and c 
are the cell edge vectors, the matrix H = {a, b,c}, and 
T is the transpose. It is obvious that using either s, or r; 
will produce the same results in physics. It further uses 
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as the particles' kinetic energy, while the true kinetic 
energy for particles is 
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where m, is the ith particle's mass. The PRMD starts 
from the Lagrangian L introduced as 


K' - $ + L 
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where <!> is the potential energy among particles and Lh 
is the Lagrangian associated with H, H, and external 
pressure and Lh has no contribution to the particles' 
dynamical equations. The aim of this comment is to 
check whether the chosen kinetic energy is appropriate 
by comparing the dynamical equations for particles with 
the Newton's Second Law 
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Case 2, let us use S{ as generalized coordinates and 
still use the complete kinetic energy of Eq. J3J). Com- 
bining Eq. ||TJ and ©, we obtain = H T F,. By 

also using of ^f-LT — miil T ^Hs, + Hs^ and -£r-K = 
m,iH T I Hs, + Hs.i ) , the Lagrangian dynamical equa- 


tion, ^ij^-L = ^-L, becomes 


m.,H J Hs, + 2Hs, + Hs, = 


H T F, 


(7) 


which is equivalent to Eq. (JSJl. 

Case 3, as in the PRMD, let us use s, as general- 
ized coordinates and use the incomplete kinetic energy 
of Eq. @. Similar to the above case, with J-if = 
TOiH T Hs'i, the particle's dynamical equation becomes 


H J Hs) 
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H T Fj, 
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which is the same as the Eq. (2) of the PRMD's orig- 
inal paper U, but different from Eq. JSJ. The dis- 
crepancy can be found in between (K — K') and 

^7 (K — K'). So the incomplete kinetic energy of Eq. (j2J 
is in-appropriate. 

If the complete kinetic energy of Eq. J2J is used, the 
first term on the right-hand side of Eq. (4) of the PRMD's 
original paperQ will disappear in physics. Anyhow, the 
PRMD is correct for crystal structure when all accelera- 
tions and velocities are zero. 
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